Enhancement of the electric dipole moment of the electron in the YbF molecule* 
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We calculate an effective efectric field on the unpaired electron in the YbF molecule. This field 
determines sensitivity of the molecular experiment to the electric dipole moment of the electron. We 
use experimental value of the spin-doubling constant 7 to estimate the admixture of the configuration 
with the hole in the 4f-shell of Ytterbium to the ground state of the molecule. This admixture reduces 
the field by 7%. Our value for the effictive field is 5.1 a.u. = 2.5 • 10^° V/cm. 



Introduction 

It is well known that effects caused by an electric dipole moment (EDM) of the electron dc are strongly enhanced 
in heavy diatomic radicals (see, for example, a review An experimental search for the EDM of the electron 

is now underway in Brighton University on the YbF molecule Q. For this reason the reliable calculations for this 
molecule are necessary. The existence of a shallow 4f-shell adds complexity to such calculations. Here we extend the 
semiempirical approach suggested in to account for a possible admixture of the configuration with the hole in 
the 4f-shell of Ytterbium. The idea that f-hole can explain small value of the spin-doubling constant 7 belongs to 
1. B. Khriplovich [|. 

The EDM experiments with radicals are made on the spin-rotational levels of the electronic ground state. The 
spin-rotational degrees of freedom of the molecule are described by the following spin-rotational Hamiltonian 



Hsr = BN^ + 7SN + SAI + WddcSn. (1) 

In this expression N is the rotational angular momentum, B is the rotational constant, S is the spin of the electron and 
I is the spin of the Yb nucleus, n is the unit vector directed along the molecular axis from Yb to F. The spin-doubling 
constant 7 characterizes the spin-rotational interaction. The axial tensor A describes magnetic hyperfine structure. 
It can be determined by two parameters: A = {A\\ + 2A±)/3 and A^ = {A\\ — A_l)/3. The last term in (|l]) correspond 
to the interaction of the EDM of the electron de with the molecular field, ■^Wd being the effective electric field on the 
electron. 

Parameters B, 7, A and are known from the experiment p0|-ff2l: 



B =. 722,1 MHz, 7 = 13 MHz, A = 7617 MHz, Ad = 102 MHz, (2) 

while Wd have to be calculated. There are three calculations of this parameter [| 13 ij]. In earlier calculations [^|jl^] 
4f-shell of Ytterbium was frozen. In jl^ f-electrons are included in the valence space, but details of this calculation 
are not yet published. Here we allow the admixture of the f-hole to the ground state of the molecule. 

Electronic wave function 

The ground state of YbF molecule is known to be S1/2 The large hyperfine constants A and Ad indicate that 
unpaired electron occupies molecular orbital CTs with dominant contribution from 6s-orbital and significant contribution 
of 6po-orbital of Yb ion. But a simple onc-configurational wave function can not explain an exceptionally small value 
of the constant 7. Following Khriplovich we will assume that there is small admixture of the f-hole: 



\T.,uj) = \(Js,ll>) +5i\ai,-uj), (3) 

where oj = ±i is the projection of the total electronic angular momentum on the molecular axis and the bar over the 
orbital corresponds to the hole. More explicitly (0) means that many electron wave function reads 
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= [...]{\af,-^)\au^)\as,Lo) +5f\ai,u;)\a,,-^)\as,^)) , (4) 

where [. . .] denotes the closed core. The spin-orbit interaction for the f-hole is large and has the opposite sign. This 
can explain the small value of 7 (see below). 

All constants of the Hamiltonian (|l|) except for the constant B depend only on the electron spin density in the 
vicinity of the heavy nucleus. For this reason we can expand molecular orbitals in spherical waves with the origin at 
the Yb nucleus: 

\as,uj) = {xs\6s) + Xp|6po) + a^dlMo) + . . .)\uj), (5) 
\af,u;) = (|4/o)+ed|Mo + ...))|c^), (6) 

where \lo) denotes spin function and dots stand for the higher spherical waves. The radial functions in this expansion 
can be considered as distorted orbitals of the Yb^ ion. Corresponding principle quantum numbers are marked with 
tilde. 

Wave function (||) is written in a pure nonrelativistic coupling case. Spin-orbit interaction Hgo mixes state (|^) 
with ni/2 states. One can expect that the largest admixtures correspond to the molecular orbitals TTp and TTf (the 
spin-orbit interaction for the 5d-shell is much smaller): 

\TTp,uj) = {ap\6p2u;) + (7) 
K,^) = (|4/2 J + ...)! -c.). (8) 

Then, the wave function of the ground state has the form 



= |cts,w) -t- CplTTp, -I- S{\a{, -u) + Cf|7rf, -uj). (9) 



Hyperfine tensor and parameter Wd 

The operator of the hyperfine interaction in atomic units has the form 



^M = |^(axr.I)l, (10) 

where g-n is the nuclear factor, a is the fine structure constant, m-p is the proton mass and a is the vector of Dirac 
matrices. This operator is known to be almost diagonal in quantum number / and in the following calculations we 
neglect nondiagonal terms. 

Let us start with the nonrelativistic expressions for a spherical wave Z 7^ 0. It is not difficult to derive, that for the 
CT-type state |/, 0) 

■^^-|g (2,'"l)(2?H-3) '"'"?l"-"^ 

So, in the nonrelativistic approximation only s-wave contributes to the isotropic constant A. Note that the radial 
integral in (^l|) is rapidly decreasing with I, while the coefficient in front of it is a weak function of I. 

Interaction of the EDM of the electron with the molecular electric field — Vt/) is also singular at the nucleus [ p5[ . 
The most convenient form of this operator is |l6j : 

Hd^2d,[ll^ (-V0). (12) 

Relativistic expressions for A, A^ and Wd can be found in |,|jl). In the paper [Q these constants were calculated 
for the wave function (^ with Cp — 5i — Cf — 0: 

A = (llSSOxf - 1024 - 54 + • • •) MHz, (13) 
Ad = (422xp + mxl + ■■■) MHz, (14) 
Wd = (-29.7xsXp + l.SxpXd H ) a.u. (15) 
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If we neglect the d-wave terms and use the experimental values (||) for A and Ad, we can solve (|l|) - (HI) for 
and Wd'- 

X, = 0.803, Xp = 0.492, (16) 
Wd = -11.7 a.u. = -1.45 • lO^^Rz/ie cm), (17) 

where sin gs o f Xs and Xp are chosen to account for the repulsion of the unpaired electron from the F~ ion. 

Res ult (p7|) is more accurate than one can expect from (^ - (^ . It follows from the proportionality between Wd 
and y/AAd'. 

Wd = -1.69 • lO^V^^d l/(e cm), (18) 

which is based on the behaviour of s- and p-wavcs in the vicinity of the nucleus where molecular potential is close to 
that of the nucleus. In contrast to that, results (|l^) depend on the assumption that spherical waves in (||) correspond 
to the orbitals of Yb+ . 

Below we calculate several corrections to ( p7| ) and (18), which mostly account for contribution of the higher spherical 
waves. If we use the normalization condition + + a;^ w 1, we can solve equations ([l^ ) - (jisl ) with the d-wave 
included: 

a;s = 0.803; Xp = 0.481, Xd = 0.356, Wd = -11-2 a.u. (19) 

Comparison of ( p7| ) and ( p^ ) gives the following correction coefficient for the d-wave contribution to Wd'. 

fed = 0.96. (20) 

Note that the normalization condition give the upper bound for the coefficient Xd. So, it is likely, that ( |l9| ) and ([20| ) 
somewhat overestimate the d-wave contribution. 



Spin-doubling constant 

Let us estimate the mixing coefficient Cp in (^): 

_ {Il{Ai)\Hso\^{X)) _ (Trp,uj\Hso\as,uj) Xpflp^e,! 



(21) 



P Ex-Ea Ex-Ea V2{Ex-EAy 

where we use (^ and (^ to calculate the numerator: 

(TTp, w|i7so|a-s, w) « Xpap{6pi,-^\Hso\6po, i) « a::pap(6pi, -i|i?so|6po, 5), (22) 

(n,Z,l,-i|i?,o|n,/,0,i) = -Cnj(',l|^i|^,0)(-i|.s_i|i) = i6uV^a + l), (23) 

and S,n,i is the atomic spin-orbit constant for the {n, Z)-shell. It is proportional to the radial integral which enters ([ill). 
For the Yb+ ion $6,1 = 1900 cm"!. 

The level (0) is identified as Ai (£'^1=18090 cm^^). It's fine sphtting with the level A2 (Hg/a) is 1370 cm^^ [0. 
Within the same approximation we can link this splitting to the constant ^g^i and find parameter Up: 

Ea, - Ea, al^e,i, ^ « 0.72. (24) 

In a same manner we can calculate coefficient Cf in (|). In this case the spin-orbit interaction mixes cxf-hole with 
TTf-hole: 

(n(j^)|g,o|I](X)) ^ (7rf,-co|g,o|af,-co) -V361U3 
Cf — ~ Of « — . 25 

Ex- Ef Ex ~ Ef Ex - Ef 

The only problem here is that molecular state with the TTf-hole, which we define as n(_F), is not known, and thus 
we do not know the energy denominator. The spin-orbit constant for the f-hole has the opposite to normal sign: 
^4 3 = —3665 cm~^. 

Formulae (^ij), ( p^ and ( p5| ) reduce the number of independent parameters in the wave function (p|) . Still, it has 
two extra parameters i5f and Sd as compared to the wave function (^ which was used in ([l^) - (p^. Below we 
eliminate parameter S{ using experimental value of the spin-doubling constant 7. 
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The spin-doubling term in the effective Hamiltonian arise from the spin-orbit interaction. It was shown in |^ 
that 

7 - 25 (1 - (EW, i| Je,+ |S(X), -i)) , (26) 
where Je = L -I- S is the total angular momentum of the electrons. For the pure E state (Je) = (S) and (|26[) gives 



7 = 0. This is no longer true, when the spin-orbit corrections are taken into account. With the help of (26), ( |2l| ) and 
( p5| ) it is easy to calculate 7 for the state (H): 

7 = 2B (2{x,a,Y-^M^ + l2Sf-MjA ■ (27) 

Note that numerical factors in parentheses are equal to 1(1 + 1). 

As we already pointed out, the second denominator in ( |27| ) is unknown. If we use corresponding energy interval for 
the Yb+ ion, we receive the following relation between Xp and 5f 

2100a;p - 11150(5f = 13, « 0.19x^. (28) 

This equation shows that experimental value of 7 correspond to almost complete cancelation between p-wave and 
f-wave contributions. To obtain this relation between 5f and Xp we used ionic spin-orbit constants S,n,i and ionic 
denominator in (^7|). So, we can not expect it to be much better than an order of magnitude estimate. 



Spin-orbit and f-hole corrections to Wd 

It follows from (|l^), (|2^) and (|2^), that for the molecular state (^) the f-wave contribution to the constant 
should be approximately 6 times smaller than that of the p-wavc. Indeed, the two contributions to the constant 7 
cancel each other, but in ( p7|) there is the factor 1(1 + 1) which is not present in (|ll|). 

A straightforward relativistic calculation result in the following f-hole correction to the hyperfine tensor 

5na ^2/^12^ 32^ 64^ \ 

" "2^ f (^i^''5/2.5/2 + ^^^^5/2.7/2 + ^^^7/2,7/2 j a.U., (29) 

fntt ,2 48 , 8 , 64 , \ 

= 2^p^' (, 245^^/^'^/^ " 49^^/^-^/^ + 441 ^^/^'^/^ ) ^'^^ ^^^^ 

The radial integrals hjj' here have the form hjji — dr {fjgj' + 9jfj'), where fj and gj are the upper and lower 
components of the radial Dirac wave function. Again we take radial integrals for the Yb+ ion and use (^8|) to arrive 
at 

SA = -llxl MHz, Md = 64x2 y^^^ j-g-^^j 



It is clear, that correction (|3l|) to isotropic constant A (|13|) is negligible, while correction to the dipole constant A^ 
(|l|) is about 15% for a given Xp, which is in a good agreement with our nonrelativistic estimate g . 
With the f-hole correction (^ included, equation ( p^ ) is changed to 

Ad = (486a;2 + 36x^) MHz, (32) 

and we obtain the following f-hole correction coefficient to Wd 

kf = 0.93. (33) 



In deriving (33) we have ignored the direct f-wave contribution to the constant Wd (see ([ISD). Such contribution is 
proportional to a small product Xd6f < 0.1 and is less than 1%. 

The last correction to the constant Wd is associated with the admixture of the TTp state to (^. According to (|2l| ) 
Cp « — 0.06a;p. This admixture changes weights of the pi/2- and p3/2-waves in the wave function. That, in turn, 
slightly changes coefficients in (|^ and (^5|) (note, that only pi/2-wave contributes to Wd)- A simple calculation give 

ho = 0.98. (34) 
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Discussion 



Taking into account (|20|), ( |33| ) and (p4[), we obtain the total correction factor for Wd to be fccorr = kdkfkso = 0.87. 
Applying it to (Jlq), we arrive at 



Wd = -10.2 a.u. = -1.26 • 10^^Hz/(e cm). (35) 

We have calculated corrections to the semiempirical value of the constant Wd of the spin-rotational Hamiltonian (|l|) 
which can account for the P, T-odd effects in YbF molecule. The main correction is caused by the f-hole admixture 
to the ground state. Two other corrections account for the spin-orbit interaction and for the d-wave term in the wave 
function. Altogether they reduce the answer by more than 10%. 

We think that wave function includes two most important configurations. Admixture of other configurations 
will tend to decrease constant Wd, so we can expect that our value (|35|) is slightly overestimated. But it is unlikely, 
that corresponding corrections are significantly larger than those discussed above. So, we estimate the accuracy of 
our calculation to be about 20%, which is typical to the semiempirical method. 

Ab initio calculation gave smaller absolute value for Wd, but it also underestimated both A and A^- It is more 
informative to compare the coefficient in the equation ( |l^ ) which correspond to different calculations. Our final value 
for this coefficient is —1.47 • 10^^ l/(e cm), while results of ||l^ correspond to —1.75 • 10^^ l/(e cm). So, in this sense, 
our result is even smaller than that of [ [l3[ . 

In this paper we have not considered other possible sources of P, T-violation, such as scalar neutral currents 0| 



and magnetic quadrupole moment of the nucleus |17|. Corresponding constants Wg and Wm of the spin-rotational 
Hamiltonian can be calculated in a similar way to Wd 0. Our value for the former constant is: Ws = —43 kHz. 
Because of the higher multipolarity of the electronic operator associated with the constant Wm, this interaction is 
much more sensitive to the higher terms of the spherical wave expansion of the molecular wave function. Thus, it is 
more difficult to make accurate calculation of this interaction within the semiempirical approach. 
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